The Loschmidt Echo as a robust decoherence quantifier for many-body systems 
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We employ the Loschmidt Echo, i.e. the signal recovered after the reversal of an evolution, to 
identify and quantify the processes contributing to decoherence. This procedure, which has been 
extensively used in single particle physics, is here employed in a spin ladder. The isolated chains 
have 1/2 spins with XY interaction and their excitations would sustain a one-body like propagation. 
One of them constitutes the controlled system S whose reversible dynamics is degraded by the 
weak coupling with the uncontrolled second chain, i.e. the environment E. The perturbative S£ 
coupling is swept through arbitrary combinations of XY and Ising like interactions, that contain 
the standard Heisenberg and dipolar ones. Different time regimes are identified for the Loschmidt 
Echo dynamics in this perturbative configuration. In particular, the exponential decay scales as 
a Fermi golden rule, where the contributions of the different S£ terms are individually evaluated 
and analyzed. Comparisons with previous analytical and numerical evaluations of decoherence 
based on the attenuation of specific interferences, show that the Loschmidt Echo is an advantageous 
decoherence quantifier at any time, regardless of the S internal dynamics. 



I. INTRODUCTION 

The physical realization of Quantum Information Pro- 
cessing (QIP) [1], requires a precise control of quantmn 
dynamics. The coherent manipulation of many-body sys- 
tems plays a crucial role for several QlP-related imple- 
mentations, such as spintronic devices 0, optical lat- 
tices 0, superconducting circuits fs'] and nitrogen- 
vacancy centers in diamond 6]. Experiments with spin 
arrays in Nuclear Magnetic Resonance (NMR) > have 
shown that many-body dynamics conspires against quan- 
tum control. Moreover, once the system interacts with 
an environment, control becomes even more difficult due 
to information leakage. This degradation of the sys- 
tem's coherent dynamics, called decoherence, is subject 
of deep theoretical and experimental investigation [lOj, 
as it remains the key obstacle for QIP. Indeed, quantum 
error correction protocols [ill, [13 can restore quantum 
information provided that they operate above a certain 
threshold. Achieving this limit is often a task for dynam- 
ical decoupling techniques [l3l - [l5| . However, specific im- 
plementations require a precise knowledge of the nature 
of the decoherence processes [l^]. It is the purpose of 
this paper to contribute to a better characterization of 
the role of a spin environment on the decoherence of a 
many-spin system. 

At least for short distance communications, spin chains 
can be used to transfer information [l3|. In fact, sev- 
eral selective polarization techniques have been devel- 
oped in NMR experiments to set up an initial local ex- 
citation in one edge of a spin chain and transfer it to 
the other edge by means of an effective XY Hamiltonian 
(i.e. SfS~+S^Sf or polarization conserving "flip-flop" 
processes) [l^, [l^. Additionally, Multiple Quantum Co- 
herence spectroscopy has allowed the study of quasi-one- 
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dimensional spin systems under the influence of spin en- 
vironments |2Cl l2l|. Here, the Double Quantum Hamil- 
tonian (i.e. the Sf S'^ +S~ S~ processes), can be mapped 
to 'AW XY Hamiltonian allowing the design and control of 
the excitation transfer in a broader family of solid-state 
spin structures ^19.. .22-25]. Thus, a deep knowledge of 
decoherence in such 1-D systems is crucial to improve the 
degree of control available for NMR-based state transfer 
protocols 26, 27]. 

A natural way to quantify the decoherence time is 
through the degradation of interferences. This requires 
the identification of specific coherence "witnesses" , such 
as excitations in the local polarization. Particularly use- 
ful are the reflections in the boundaries that can be ob- 
served as well defined Mesoscopic Echoes (ME) (28l - [30| . 
Recently, the ME intensity has been used to quantify 
decoherence of spins arranged in a ladder topology |31| . 
Alternatively, the evaluation of can be performed by a 
time reversal procedure, the Loschmidt Echo (LE) [33], 
where one evaluates the reversibility of the system's dy- 
namics in the presence of an uncontrolled environment. 
The LE can be accessed experimentally in many situa- 
tions, such as spin systems [3i[8ll33|, confined atoms [13] 
and microwave excitations (35|- Besides, it has become 
a standard way to quantify decoherence, stability and 
complexity in dynamical processes, in several physical 
situations [s^, [33] . 

In the present article, we address the controllability of 
a spin chain (5) in the presence of a spin environment 
{£) by performing a quantitative study of the LE. The 
LE degradation characterizes the decoherence due to the 
perturbation of £ on the otherwise simple dynamics of 
S. Indeed, the many-body nature of the S-£ interaction 
yields a very rich behavior in the dynamical regimes of 
decoherence: a short time quadratic decay, an exponen- 
tial regime, and a saturation plateau are identified by 
our numerical approach. In particular, we perform a de- 
tailed analysis of the LE exponential decay, addressing 
how the rates scale with a Fermi golden rule (PGR). Ad- 
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ditionally, since for weak perturbations the LE of the lo- 
cal excitation can be seen as a Survival Probability (SP), 
the numerical results are compared to previous analytical 
predictions for that magnitude 

In the next section, we describe the spin problem and 
summarize its theoretical background. We introduce the 
spin autocorrelation function and describe a single parti- 
cle analogy which underlies further analysis. We also dis- 
cuss the local version of the Loschmidt Echo and present 
its definition in terms of the local polarization, which 
is the usual experimental observable. In Section IIII A| 
we present the numerical study of the LE for some of 
the different physically relevant parameters. We con- 
sider S-£ interactions which are weak compared with 
those determining the £ dynamics. We also address dif- 
ferent anisotropics of the S-£ interaction which ranges 
from pure XY (planar) to truncated dipolar cases, going 
through the Heisenberg (isotropic) interaction. In Sec- 
tion |TiTB] we analyze the obtained results. First we con- 
sider the transition from the short time quadratic decay 
to the exponential regime in analogy to what is known for 
the SP. Then, we focus on the exponential regime to show 
that the Fermi golden rule in the present spin problem 
has independent contributions arising from each specific 
process in the S-£ interaction. We then compare these 
rates with a previous evaluation based on the contrast 
of specific interferences (ME attenuation). In the last 
section, we conclude that since the LE filters the internal 
dynamics of the system, it provides a reliable and contin- 
uous access to all time regimes. Thus, the LE compares 
favorably with the evaluation of decoherence based on 
interference contrast. 




FIG. 1. (Color Online). The spin system, (a) Open boundary 
conditions, (b) Closed boundary conditions (ring-like). Con- 
tinuous (green) connections represent interactions that can 
be inverted to obtain the Loschmidt echo. Dash (blue) lines 
represent non-controllable interactions. The first spin (black 
circle) is initially polarized and the rest of the spins are in one 
of the high temperature configurations. 

flip-flop interaction between nearest neighbor spins. In 
the model of Fig. [l]-a, i.e. the chain: 



II. QUANTUM DYNAMICS IN SPIN SYSTEMS 
A. THE MODELS 

The spin models analyzed in this article are schema- 
tized in Fig. [TJ In the first one, the system S is an m-spin 
chain (Fig. [l]-a), which could constitute a quantum chan- 
nel. It interacts with a second chain £, that stands for the 
"environment" which perturbs the dynamics of S. The 
second model (Fig. [T]-b) is obtained from the first one by 
imposing a periodic boundary condition that transforms 
chains into rings. 

For both models, the spin Hamiltonian is given by: 

Hiatal ^Hs(^h+IS(^H£ + Vs£ , (1) 

where the first and second terms represent the system 
and the environment Hamiltonians respectively, and the 
third one is the interaction between them. In order to 
simplify notation, we write just and instead of the 
tensor product with their respective Ig and Ig identities. 
For both i' ~ S or £, we use an effective "planar" or 
XY Hamiltonian that describes the homogeneous 



m — 1 

71=1 

m — 1 J 

— ^ ~^i^L^,n+l^iy/n ^i^,n+l^}!^,n)- (2) 

Here S^^^ and iS^^ are the x and y components of the 
spin operator in the n-th site in the v chain respec- 
tively, while and are the raising and lower- 
ing operators. Again, the abbreviated notation for any 
spin operator must be understood in the form „ = 

A ® . . . (8) 5^ (8) . • . <8) /m (E) Im+i (8i . • . <8) /2m • If ouc wants 
to consider the ring model, an extra XY coupling ap- 
pears between the 1-st and m-th spins. 
The interchain coupling is: 

m 
n=l 

(3) 

m ^ 

n=l 

(4) 



3 



where the first term is an Ising interaction. The a pa- 
rameter determines the anisotropy of the coupling. This 
encompass the typical magnetic resonance scenarios: the 
XY (planar) interaction [18j, represented by a = 0; the 
Heisenberg (isotropic) interaction 39], by a = and 
the truncated dipolar coupling corresponds to a = 1. 
In order to extend and systematize our analysis we also 
consider several other values for a. It is important to no- 
tice that for any finite a the S-£ interaction has always 
an XY component. This allows a polarization exchange 
which, in a Fermionic representation, can be seen as a 
"single-particle tunnelling" [4l|. In such a picture, the 
Ising term corresponds to a nearest neighbour Hubbard 
term which is a two-body interaction. 

It is crucial to stress that the real constants Jj, 
and Jjg determine the relevant time scales of the whole 
dynamics. As introduced above, the first two give the 
homogeneous XY coupling within S and £ respectively, 
while Jgg stands for the interchain coupling. To ensure a 
smooth degradation of the S coherent dynamics, we set 
Jjg in the weak coupling limit, i.e. J^^ <^ Jg = Jg. 



see Appendix |X] for further details. Here, each contri- 
bution I^E"^) to the locally polarized initial state can be 
written as: 



|*,) = |ti)®|/3.), (7) 
where the basis for the remaining 2to — 1 spins is 

|/3r) = \s2)®\s-i)®\si)®...®\s2m) with |sfc) G {|t) , ■ 

(8) 

Since in the regime of NMR spin dynamics, the thermal 
energy ksT is much higher than any other energy scale of 
the system d^Sl], all statistical weights result identical, 
i.e. \cr\ = 

It is useful to analyze the autocorrelation ^ (t) . For 
each contribution to the initial state, the evolved 
wave function at time t, is 



|vE'.(t))=exp{-liJ,„,,,i}|*,.) 



(9) 



Then, the probability of finding the first spin up- 
polarized is: 



B. Basic features of Spin Dynamics 

A natural question for the spin models introduced 
above is how to quantify the decoherence of S in the 
presence of £. This means that one has to deal with a 
composite (bipartite) Hilbert space S ^ £, and trace out 
the ^-degrees of freedom whenever necessary. A stan- 
dard strategy relies on finding an appropriate decoher- 
ence rate 1 /r^ through the quantification of the attenu- 
ation of system's specific interferences. As in the experi- 
ments [1^, [1^ , one starts evaluating the evolution of an 
injected local polarization through the spin autocorrela- 
tion function [33, B^l : 



Pi,i{t) = 



(^e,|%iW%i(0)|^e,) 
(*ej%i(0)%i(0)|*,,) 



(5) 



This function gives the local polarization at time t along 
the z direction in site 1 provided that at time t — 
the system was in its thermal equilibrium state plus 
a local excitation in site 1. Here, the spin operator 
in the Heisenberg representation is given by Sg ^ (t) — 



5,16"'^*°*"'*/''. The many-body state \^ ^q) 
corresponding to a high temperature thermal equilibrium 
represents a mixture of all states with amplitudes sat- 
isfying the appropriate statistical weights and random 
phases. Then, the initial non-equilibrium local excita- 
tion \^ne) can be defined in terms of the computational 
(Ising) basis components that have the I''' spin up as: 



1^, 



1/2 



(6) 



^^l:iw = Ei((tii®(/3.i)i*.w>r 



(10) 



where the sum runs over the 2^™~^ configurations of the 
2m — 1 remaining spins. Notice that the sum over the 
j-index in Eq. I10[ means that we are performing a trace 
not only over f , but also over the spins in S other than 
the first one. After the summation over all contributions 
\^r) to the initial state, and expressing the result as a 
local polarization (2^, we obtain: 



Pi At) 



E 



X 2. 



(11) 



A computation of the time dependent local observable in 
Eq. [5] reduces to Eq. Ill I which in turn, requires evolving 
each of the 2'^'"^^ pure states to evaluate the ensem- 
ble averaged observables. This is implemented using a 
Trotter decomposition [l^l assisted by an algorithm that 
exploits quantum parallelism [4^ (See appendix [BJ . 



The single-particle picture and Mesoscopic 
Echoes 



The obvious complexity of the many spin dynamics 
might hinder some simple interference phenomena that 
can be taken as hallmarks. Under certain experimentally 
achievable conditions (e.g. XY interaction, ID topol- 
ogy and negligible S-E interaction), the autocorrelation 
Pi^i(t) becomes a simple one-body magnitude. Indeed, 
once the initial excitation is created, the first physical 
picture about its evolution may be obtained from the 
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Wigner- Jordan spin-fermion mapping [46|]. The point 
here is to trace over the ^-degrees of freedom from the 
beginning and focus on a single quantum spin chain. Ac- 
cordingly, an isolated m spin chain with XY interaction, 
where N of them are up, is mapped to a chain with 
N non-interacting fermions. Thus, a local polarization 
excitation has the same dynamics as a single fermion 
wave-packet in a tight-binding linear chain [4l|, \^ . The 
observed autocorrelation function in the limit of infinite 
temperature is precisely described by the evolution of a 
single spin up in a chain of down spins: 



Iti) ® 112) ® \U) «) ... «) \u 



(12) 



This is a one-body wave function, defined on a sub- 
space of S where the total spin projection is ^ ~ 1 ; 
the observable is evaluated as: 



(13) 



Since this is a finite size system one should expect revivals 
of the initial polarization. Such revivals are called Meso- 
scopic Echoes (ME) [2^ [s^l and they appear when con- 
structive interferences manifest at the Heisenberg time 
tn ~ ^/A, with A being the typical mean energy level 
spacing. In a spin chain with XY interaction one may 
safely use A ~ Jg/m. As long as such one-body picture 
remains approximately valid for a linear chain weakly 
coupled to the environment these interferences show up 
experimentally [H, [H, ill . 

A weak coupling to the spin bath £ results in a pro- 
gressive attenuation of the ME, which has been used to 
quantify environmentally induced decoherence (3lj | . This 
attenuation is understood as a Fermi golden rule (PGR) , 
which describes an "irreversible" decay of a pure state 
in S into collective S ^ £ states. In fact, the validity of 
the FGR requires here the breakdown of degeneracies, 
i.e. the whole S ^ £ must behave as a fully many-body 
system. We shall return to this point below. 



to a non-reversed perturbation. 



(14) 



acting in both periods. Thus, in analogy with Eq. [SJ we 
define the observable Loschmidt echo as the recovered 
local polarization: 



MLE{2tR) 



(^e,|^il(2^i^.)%i(0)|^e,) 
(*e,|%i(0)%i(0)|*e,) ■ 



(15) 



The spin operators, expressed in the Heisenberg repre- 
sentation, are now: 



(16) 

The computation of the Loschmidt Echo in Eq. [15] pro- 
ceeds as we did above for the forward dynamics. In a 
system with a general interaction, it requires a full many- 
body evolution. 

As before, the initial excitation is described in terms 
of the Ising basis by Eq. [6l For each contribution |^',.) 
to the initial state, the resulting wave function at time 
t = 2tj^, after the whole time-reversal procedure, is 

|*,(2tfl)) = exp{-l(-7?5+S)tfl,}exp{-i(iJ5+S)tfl.} j*,) . 

,,(17) 

In analogy with the discussion above, the probability of 
finding the first spin up-polarized is: 

<lW=EK(ti|®(/3,l)l*.(0>l', (18) 



where the sum runs over the 2^™"^ configurations of the 
2rn — 1 remaining spins. Again, the sum over j means a 
trace over all the spins of the environment and the system 
except the first one. After the summation (average) over 
all contributions j^'^) to the initial state, and expressing 
the result as a local polarization, we compute the local 
Loschmidt Echo as: 



D. Loschmidt Echo 

Let us now explain the essence of the protocol that 
uses the Loschmidt Echo to quantify decoherence in a 
spin system using a local spin as an observable Q- Our 
strategy relies on the controllability of the chain S, whose 
Hamiltonian's sign can be switched at will as it is often 
the case in NMR. The reversibility of the dynamics within 
the chain is perturbed by the interaction with the non- 
controlled spin chain £. There are two stages in the evo- 
lution. First, a spin excitation is created and the whole 
spin set evolves according to the Hamiltonian of Eq. [TJ 
during a time tR. At that time, the internal interactions 
within S are reversed (i.e. is replaced by —Hg during 
a second i/j-period). However, neither the S-£ coupling 
nor the interactions within chain £ are reversed, leading 



X 2. 



(19) 



Here again, the statistical weights \cr\ are the inverse of 
22m-i^ the number of initial states in the ensemble that 
satisfy the "1''* spin up polarized^^ condition. Notice that, 
except for the fact that the evolution operator contains 
a partially reversed dynamics, this quantity refers to the 
same physical observable as Eq. [TTl 



E. Connection to previous works 

It would be useful to make a connection between the 
observable just described and the usual definition of the 
Loschmidt Echo [3^ . which was extensively studied in 
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the field of quantum chaos and quantum information 
[sgI [37} . With this purpose, let us consider the particu- 
lar case of Fig. [1] where describes a chain with XY 
interactions, while describes a chain £ that remains 
quenched in a random configuration and V^^ is restricted 
to an Ising interaction. Under these assumptions, the 
time-reversed dynamics reduces to that of a single spin 
up in an oriented chain defined in Eq. [1^ Accordingly, 
S becomes an Hermitian self-energy operator £5 act- 
ing on the S Hilbert space. Indeed, S5 represents a set 
of non-reversed random energy shifts produced by the 
Ising interaction with the static environmental spins. In 



MLE{2tR) = 



Here, one can recognize the LE definition introduced in 
Ref. [S^I as the overlap of two wave functions evolving in 
presence of a quenched disorder whose dynamics is not 
inverted. 

The decomposition of the spin set into a controllable 
subset {S) and an uncontrollable one {£) resembles the 
discussion of the partial fidelity, called Boltzmann echo, 
analyzed in Ref. [sOj for a two-body problem. Analo- 
gously, the 2m spin problem treated here verifies that: (i) 
a separation into two interacting subsystems [S and £) 
is performed, (ii) the initial state in S is at least partially 
prepared (injected polarization), and at the end a local 
measure is performed in the same site of injection, (iii) 
the subsystem £ remains in the high temperature ther- 
mal equilibrium, and (iv) the Hamiltonian in S can be 
fully time-reversed, while the Hamiltonian of £ and the 
S-£ interaction remain uncontrolled. However, the main 
result of Ref. [SO] , that considers a chaotic one-body sys- 
tem coupled to a chaotic one-body environment, can not 
be directly compared with our study. Here, we focus on 
a system and an environment which are both many-body 
systems that can be reduced to integrable one-body sys- 
tems. In our case, we expect many-body chaos |5l| only 
as a consequence of the S-£ interaction. 

Many-body chaos has been the subject of much interest 
[s^ - fs^ l , which was recently renewed mainly in connection 
to thermalization dynamics [ssl j. Within this context, 
the attention has been centered on the study of spectral 
correlations and the related survival probability (SP), i.e. 
the decay of an initial excitation [3g]. Based on general 
considerations on the strength function or local density 
of states (LDoS) t56|, one may predict a Gaussian (or 
quadratic) time decay governed by the second moment of 
the perturbation cr^, followed by an exponential regime 
whose rate is described by the PGR: 

-^a^iVi. (21) 



the independent fermion picture, E5 yields a "random 
potential" for each specific configuration of the environ- 
ment, i.e. the the binary alloy variant of an Anderson's 
disorder [4^. It is now relevant to address the physi- 
cal meaning of the sum over the 2^™"^ indices r of the 
observable in Eq. 1191 This sum performs a trace over 
system spins which evolve but are not observed (crucial 
to recover a one-body dynamics), as well as a trace over 
the environmental spins, which can be seen as an ensem- 
ble average over quenched disorder configurations [3^ . 
Thus, the same procedure that enabled to reduce Eq. [5] 
to Eq. [T31 transforms Eq. [15] into the corresponding 
one-body Loschmidt Echo: 



(20) 

I 

Here, Ni is the density of directly connected states 
(DDGS). While it might seem elusive, this magnitude 
can be numerically evaluated with a Lanczos algorithm 
(57| . It also accounts for the transition time from the 
short-time quadratic regime to the exponential one. Ac- 
cordingly, the cross-over to the exponential regime is ex- 
pected [allli] to occur at ts ~ [TpGR'^'^y^ = {2TT/h)Ni. 
This spreading time characterizes the dissemination of 
an excitation within the environment. In the many-body 
context it has been proposed an interpolation formula 

Hi, 

Pn{t) = exp[2^ - 2^^ - THyh% (22) 

which, to our knowledge, has not been checked in con- 
crete systems. 

Goming back to our model, when the chain S is iso- 
lated, the localized excitation propagates freely, and the 
LE would have a constant value of 1. However, once a 
weak Vgg interaction is turned on, the excitation decays 
into the chain £ and the local LE degrades progressively 
with a law that would be closely described as a SP. Thus, 
for long chains where the spectrum is a quasi-continuum, 
we expect that the short and intermediate time-regimes 
of the LE would follow closely the above discussion for 
the SP. In fact, we analyze the LE decay numerically and 
study the rates testing the validity of a PGR description 
[59I]. When applied to the present many-body context it 
would look: 

where jV^^I^ = is the local second moment of the 
process S (e.g. Ising or XY) contributing to the S-£ 
interaction and Nis represents some appropriate DDGS. 



exp{-^(-i?5 + ±s)tR}cxp{~^{Hs + ±s)tR} 1*1) 
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[units of hli^ 



FIG. 2. (Color Online). Time regimes of the local Loschmidt 
echo. Numerical results for a ring of 5 spins weakly cou- 
pled to another identical ring by an XY interchain interac- 
tion (a = 0). The Mle is plotted as a function of the total 
evolution time t — 2tR. Notice that the stronger the inter- 
chain coupling Jgg , the faster the saturation regime is reached 
(dashed line). From up to down, the different curves corre- 
spond to the coupling parameters Jgg: 0.001, 0.01, 0.025, 
0.05, 0.1, 0.25 and 0.5, in units of Jg. 



III. QUANTIFYING DECOHERENCE 
THROUGH THE LOSCHMIDT ECHO 

A. Numerical Results in Spin Chains 

In this section we present the results obtained for 
MLE{t) in the spin models represented in Fig. [T] Even 
though our major concern lies on the exponential decay 
described by the FGR, we can also identify the short time 
quadratic decay and the saturation regime, as shown in 
Fig. [2l It is noticeable that the Loschmidt echo yields 
results for a very wide range of parameters and times 
scales. This feature contrasts with the study of interfer- 
ences through the ME whose observability restricts the 
quantification of decoherence. 

In Fig. [3]-a the short time dynamics is compared with 
the expected quadratic decay (dotted line), which should 
appear in any quantum evolution at short times. Indeed, 
the plot of (1 — Mle) / Jse (^'^ log- log scale) as a function 
of time, shows that Mle follows a quadratic function: 



MLE{t) ~ 1 



(24) 



This confirms that the short time decay is determined by 
CT^ = (Jgg/2H)^, the second moment of the S£ interac- 
tion. The agreement is observed until a time tg ~ h/ Jg, 
which verifies the prediction for ts in terms of the £ dy- 
namics. For comparison, we show in Fig. |31-b the SP of 
an excitation in a single spin system S that interacts with 
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FIG. 3. (Color online). a)Short time behavior of the 
Loschmidt Echo. The dotted line is the quadratic decay (Eq. 
1241) . b) Short time behavior of the Survival Probability for an 

XY spin chain described by: H total ~ Jse (^^1 ^2 + ^2 ^1^ + 

), where the first spin is the sys- 
tem and the others constitute the environment. This system is 
known to satisfy the FGR 6Q| . Notice the similarity between 
(a) and (b) in departing from the quadratic decay, which 
evinces the spreading time t^. The dashed line shows a plot 
of Eq. [22] for Jqe = 0.5 Jg. 



the edge of a spin chain £ where it can spread through 
pure XY interactions. This model, shown in the inset, 
constitutes a paradigm for the onset of the FGR, because 
the DDCS is precisely defined by A^i = 1/ Jf and remains 
independent of Jg^ [60|. In Fig. |3]-b, we also show with a 
dashed line the interpolative expression in Eq. [2H for the 
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20 30 40 50 
t [units of h/Jn] 



FIG. 4. (Color online). The onset of the LE exponential de- 
cay. In the log scale, the dashed lines are the ideal exponen- 
tials determined by the fitting, show here for representative 
values. The delay in the entrance into the exponential regime 
is indicated with arrows pointing the first data used in the 
fitting. Continuous lines represent the interpolative Eq. 1221 
used with the corresponding and F. 



strongest coupling. This last expression deviates faster 
from the quadratic decay than the SP in the one-body 
chain. 

The onset of the exponential regime is exemplified in 
Fig. m for a few values of a pure XY S£ interaction. 
As a general tendency, we observe that the asymptotic 
exponential decay only becomes well defined after some 
evolution period. We indicate with an arrow the initial 
data that we use to fit these rates. For comparison, the 
SP interpolative expression is plotted for the same decay 
rates. We notice that for Jgg/ Jg ^ 0.085 the interpola- 
tion starts to lie below the numerical results, and longer 
times are required to define the asymptotic exponential. 

For long times, the LE shows a saturation plateau at 
1/10 (see Fig. [5]). Such observation is consistent with the 
expectation that, within these coupling networks, the fi- 
nite system of interacting spins behaves ergodically under 
the Loschmidt Echo dynamics, and thus the polarization 
spreads uniformly. At long times each site is polarized 
by the amount l/(2m)The larger the coupling J^^, the 
faster the asymptotic saturation is reached. 

In order to assess the exponential regime, we plot the 
characteristic rate l/r^ in Fig. [5]-a, as a function of Jjg 
in units of J^/h. This quantity is appropriate to verify 



the FGR (Eq. [23]), as long as Jjg is the typical scale for 
the second moment of the S-£ interaction and 1/ that 
of A'^i, the DDCS. We point out that rate calculations 
with longer rings and chains are seen to give similar val- 
ues as long as the above parameters are kept the same. 
Such independence on m resembles to the SP and LE for 
a single spin interacting with chains of different lengths 
[goI [61| . This is indicative that we are reaching the limit 



were the DDCS becomes dense enough to manifest the 
sum rule intrinsic to the FGR[5^. Thus, we restrict our 
analysis to cases with m — h. Although several forms of 
the interchain Hamiltonians were considered by varying 
the anisotropy parameter a , we show only those relevant 
to NMR experiments: XY {a = 0), isotropic {a = —5), 
and truncated dipolar interaction (a = 1). We observe 
that the boundary conditions play a non-trivial role f4l[ . 
For the case of open boundary conditions (Fig. [T]-a), 
some oscillations appear mounted on a decay which also 
depend on the parity of m, the number of spins in each 
chain. Here, we present only the results using closed 
boundary conditions (rings) where these effects are much 
weaker. 

Note that the Loschmidt Echo allows one to explore 
a range of very weak perturbations yielding from the 
regime of exponential decay. Such a ran ge w as inaccessi- 
ble in the interference contrast method [3l|. We observe 
that the rates in Fig. [5]-a start from zero and have a 
rapid increase with the Jse perturbation. After a small 
perturbation threshold the rates slow down to a linear 
dependence on the second moment of the perturbation. 
This confirms the validity of Eq. [53] for this range. The 
linear fit is shifted by an offset, I/tq, which seems to de- 
pend on the nature of the S£ interaction as it becomes 
larger for those perturbations with bigger Ising compo- 
nents. 

Fig. [5]-b shows the FGR contribution to the scaled 
decoherence rates [I/t^ — l/To]/Jgg as a function of in- 
teraction anisotropy a^. There, we also include the rates 
derived from the ME attenuation [3l| ■ From the slope in 
that plot we derive the contributions to the global decay 
rate I/t^, arising from XY and Ising processes in the 
interchain interaction; 



1 



1 

To 



'4> 



For the Loschmidt Echo results: 
1 



LE 



LE 



72 

(0.92 ±0.04)-^, 



(25) 



(26) 



-^ = (1.12±0.04)a2^. (27) 



In order to compare with the numerical results in Ref. 
[31], and translating its notation —a/b = 2a, the rates 
contributing to the ME degradation result: 



1 P 
ME: = (1.00 ±0.06)^, 



(28) 



1 P 
ME: — (2.0±0.3)a^-^. (29) 



^zz 
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FIG. 5. (Color Online), a) Decay rates for the NMR relevant interchain couplings, in units of Jg/h. The slopes and offsets 
depend on the value of ct, i.e. on the relative weight between the Ising (dephasing) and XY (polarization transfer) contributions 
to the interchain coupling, b) Scaled rates showing the additivity of the XY and Ising contributions to the FGR rate as a 
function of . For comparison, the rates obtained from Mesoscopic Echo degradation data obtained in Ref. [sJl are also 
plotted. 



The most striking effect evinced by this comparison is 
that while the XY contributions are essentially the same 
in both methods, the Ising contribution (i.e. pure de- 
phasing) in the ME (Eq. is almost twice the value 
obtained from the local LE. In the following section we 
analyze the origin of such result. 



B. Analysis of the Decoherence Regimes 

First, we would like to address the applicability of a 
"iSP picture" for the LE in the present case, as it was 
proposed in Section|TT]in connection with previous works. 
We verified that the short time decay is naturally ruled by 
the second moment of the S-£ interaction = 
We also checked the estimation for the spreading time 
ts — fi/Jg, *.e- the transition from the quadratic short- 
time to the exponential regime. We found that such 
prediction is in agreement with the observed end of the 
quadratic decay. 

In our spin problem, the LE exponential regime needs 
a further delay beyond to become well defined. This 
indicates that the DDCS, is not immediately defined but 
Jjg may have a role in its stabilization. Indeed, it is 
needed a back action of J^g to break the strong degen- 
eracy present in the unperturbed integrable f , a finite 
XY chain itself. Only then, one has the sufficiently dense 
spectrum required for the FGR to apply. Once the whole 
Hilbert space is made accessible by the perturbation, as 
many as (^^) states become coupled by an interaction 
that conserves spin projection. 

The interpolative Eq. [51] approximates satisfactorily 
the short-time decay, but it is not as well suited for in- 
termediate times and certainly can not be applied to the 
saturation regime. These deviations can be, in principle, 
related to a failure in the SP picture to describe the local 



LE dynamics. The physics underlying Eq. [H] is based 
on the assumption of a sufficiently large number of par- 
ticles and a well defined A^i. Also, one should be aware 
that our spin case is strictly finite, indeed quite small, 
and therefore the polarization cannot relax to zero. The 
observed asymptotic steady magnetization 1/2™ can be 
identified with an ergodic behavior of the excitation de- 
scribed by the local LE dynamics. Since this ergodicity is 
not present in the unperturbed system, it should emerge 
as a consequence of the S£ interaction. 

It is important to notice that neither the rate obtained 
nor depend on the total number of spins or on the 
number of spins in S. The reason for such independence 
is related to the initial non-equilibrium condition stated 
in Eq. [51 which is a well localized excitation that main- 
tains its character when it propagates through an XY 
chain. It turns out that such localized excitation behaves 
much like a single particle and, in some sense, almost clas- 
sically. This remains true at least when the LE is in the 
exponential regime. We have seen that the independence 
on the number of spins does not hold if one considers an 
initial state built-in as an arbitrary superposition. How- 
ever, the investigation of this issue goes beyond the scope 
of the present study. 

One should notice that splitting the FGR contribution 
to the decoherence rate into two well separated terms, a 
strategy also exploited in Ref. (sij, evidences the addi- 
tivity of the XY and Ising processes. Each of them is 
associated with a different term in the interaction cou- 
pling and require a different DDCS. These properties 
are indeed manifested in Fig. [5l that shows the linear de- 
pendence of (1/t0 — l/To)/Jg£ on the relative weight a^. 

Comparing the rates obtained through the LE with 
those computed by the ME degradation 31], we notice 
that the XY rates are equivalent and we can interpret 
such equivalence by using the mapping into a one-body 



evolution in both S and £. Even when the one-body 
picture is not rigorously valid when the interaction Vgg 
is turned on, it provides some physical insights (6lj | that 
apply to more complex cases. Accordingly, the dynamics 
along the chains is only weakly affected by the tunnelling 
processes, (i.e. in a single particle picture, the kinetic 
energy along the chains commutes with that along the 
interchain direction). Thus, the rate should co- 

incide with the interchain tunnelling rate and we expect 
that it should not be affected by a time-reversal proce- 
dure within S. 

In contrast to the simple decay associated to the XY 
component of the S-£ interaction, the Ising part causes 
energy fluctuations that induce a diffusion-like process 
within the system. This tends to blur out the dynami- 
cal recurrences (i.e. MEs). The smaller decoherence rate 
observed from the LE indicates that the time reversal 
procedure is at least partially effective to reverse such 
processes. In fact, the rates for the pure Ising interac- 
tion observed here from the LE are about half of those 
obtained from the ME [3l|. This means that the ME at- 
tenuation overestimates the phase degradation induced 
by the environment. As stated before, the displacement 
of a spin excitation in presence of a quenched spin envi- 
ronment plays the role of an Anderson disorder degrading 
the wave packet dynamics that produces the ME. Such 
degradation is computed as decoherence. The reversal 
of the internal XY interaction would not be able to re- 
verse such disorder. However, if the environment has 
an inner dynamics in a time scale comparable to that of 
the system {Jg ~ J^), there are particular fluctuations 
that allow for a perfect reversal of the system dynamics. 
Indeed, this occurs when both the hopping and the site 
energy signs are inverted. Such specific ^-fluctuations are 
those needed to unravel the phase shifts produced during 
the forward evolution. This argument, which relies on a 
single particle picture and equivalent time scales for the 
local dynamics in S and £, is represented in Fig. [51 Thus, 
one may safely say that in the presence of a fluctuating 
environment, the LE is able to reconstruct the phases in 
an appreciable fraction of the local configurations. 



IV. CONCLUSIONS 

We presented an evaluation of decoherence for a spin 
chain in realistic many-body scenarios. As specific re- 
alizations for structured environments, we used a sec- 
ond spin chain weakly coupled to the first. This system- 
environment interaction ranges from pure XY to trun- 
cated dipolar, passing through the isotropic Heisenberg 
interaction. In order to evaluate decoherence, we resorted 
to the Loschmidt Echo, which here is the local polariza- 
tion recovered after an imperfect time reversal procedure. 
The attenuation of such echoes yields a first hand estima- 
tion of the decoherence rate, without any ad-hoc assump- 
tions abo ut sp ectral functions [6^ . Issj or stochastic noise 
operators |64l. l65j . Our computational technique involves 
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a) Forward 

(>-iP><)5 

b) Imperfect Backward c) Perfect Backward 

I I I I till 

FIG. 6. Spin propagation in presence of an Ising interchain 
interaction, represented as a single fermion moving through 
a binary alloy landscape. A filled circle is a spin up or a 
fermion, and an empty circle represents a spin down or a 
hole. All the processes are considered in the short time scale 
given by a single Trotter time step, a) Forward dynamics 
within the upper chain <S, in presence of the "random site 
energies" produced by the Ising interaction with the lower 
chain £. b) If £ remains frozen, the backward evolution in 5 
is imperfect because the signs of the "site energies" are not 
inverted. c)A particular evolution of £ provides a perfectly 
reversed dynamics for S. In this case, the local portion of 
potential landscape, that determined the excitation's forward 
evolution, is reversed. 



a Trotter Suzuki decomposition assisted by a recently de- 
veloped algorithm that relies on quantum parallelism to 
evaluate local observables (isf . 

In the present many-body context, the decoherence 
rate is separated in two contributions, both of which scale 
with Jj^ and 1/ Je, i.e. as a Fermi golden rule. The rates 
obtained here do not depend on the number m of spins. 
This is strongly related to the localized initial condition 
but such independence does not hold if the initial con- 
dition is a superposition state. It is also indicative of 
a specific sum rule relating the local second moment of 
the interaction and the DDCS, which in this many-body 
case coincides with that resulting from a single particle 
picture. 

In the adopted model, the close connection between 
the Loschmidt Echo and the Survival Probability allows 
for a comparison with the expectancies of the last magni- 
tude. In particular, we assessed an interpolative formula 
[3^ that matches the initial quadratic decay with the ex- 
ponential regime. Despite of the qualitative agreement, 
the Loschmidt Echo evidences a richer and more complex 
dynamical behavior than those predictions. For instance, 
the breakdown of single-particle degeneracies due to the 
S-S interaction and the appearance of an ergodic regime, 
manifested as a steady saturation, are now clearly shown 
in the numerical results. 

The numerical studies performed here confirm that the 
Loschmidt Echo is better to quantify decoherence than 
the standard analysis based on interference degradation, 
as it recovers information that escaped such analysis [sij . 
Additionally, by compensating the intrinsic dynamical in- 
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terferences of the system through time reversal, the LE 
gets rid of most of the trivial part in the S dynamics that 
conceals the decoherence effects. Thus the LE provides 
a smooth and continuous access to characterize decoher- 
ence processes. 
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Appendix A: The spin autocorrelation function. 



We show here equivalent expressions for the spin autocorrelation expressed in Eq. [5] The autocorrelation Pi,i(t) 
gives the local polarization at time t along the z direction in site 1 provided that at time t = the system was in its 
thermal equilibrium state plus a local excitation in site 1. In order to show this assertion explicitly, we rewrite the 
numerator in Eq. [S] 



{^JSt{t)Sm I*,,) = (^S+it)S^it) - (^^+(0)5r(0) - I*,,) , 



(Al) 



where we have used the identity — S^S — i, which is valid for spins 1/2. The S label has been dropped for 
simplicity. Then, 

1 



(A2) 



where we have rearranged the terms and used the high temperature hypothesis for j^gg) explicitly, i.e. 
{'^eq\ ^i{t)\'^eq) — ("^ eq\ ^ 1 i^) l"^ eq) — 0- I* is crucial uow to identify the meaning of the expectation value 
{^eq\ ■■ ■ I'^eq) ^ tracc ovcr the basis states (see Eq. (3) in Ref. [4l|). Hence, the invariance of the trace un- 
der cyclic permutation leads to: 

1 



(A3) 



and one identifies the non-equilibrium initial condition S^{0) |^eg): which after normalization yields exactly Eq. [HI 



Accordingly, the numerator can be written as: 



{'^eq\SsJl^\^eq) 



1/2- 



(A4) 



= ^(*„e|5+(t)5r(t)|f,J-i. 

The denominator in Eq. [5] is (^^J Sf{0)Sf{0) j^^g) = j- Therefore, 

Pia(i)=2(vI.„J5+(t)^r(0|*„,)-l. 

In the Schrodinger picture it yields: 



Pi At) 


= 2 




iH, , 




Pi At) 


= 2 






{s+s^ 


Pi At) 


= 2 


(*„e 


iH, , 





I ^ nel I 



(A5) 

(A6) 

(A7) 
(A8) 
(A9) 
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Thus, the autocorrelation written in the form of Eq. IA9l is. in fact, the expectation vahie of the operator 

over the evolved non-equilibrium state e~'^*°*°'*/'' |^'„e)- We note that an equivalent reasoning can be performed in 
the case of the local Loschmidt Echo in Eq. [15] just replacing the forward propagator 6"'^*°*"'*/'' by the Loschmidt 
Echo propagator e"i:("-f^s+^)*«e"S(^'S+^)*«. 

Appendix B: Numerical Simulations with Trotter-Suzuki evolutions. 

The computation of the time dependent observable in Eq. [S] requires to evolve each of the states participating in the 
ensemble. However, its implementation for large systems, e.g. by means of a 4*^* order Trotter-Suzuki decomposition 
(4^ . soon becomes quite expensive in computational resources. Instead, the local nature of the excitation allows the 
employment of an algorithm exploiting the quantum parallelism [451] to give an exponential reduction of computational 
efforts. The evolution of a few pure states is enough to describe the ensemble averaged excitation dynamics. Thus, a 
typical initial state representing the whole ensemble has the form: 

|*> = Iti) ® I E ^/^jL^T ^'"" l^^^l ' = random phase, (Bl) 

which exploits the quantum superposition over the whole spin set. Typically, a few of these entangled states is enough 
to get rid of statistical noise and obtain local physical observables with good accuracy. This highly efficient technique 
for spin-ensemble calculations is enhanced by the parallelization enabled by its implementation on a General Purpose 
Graphical Processing Unit (GPGPU). 
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